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ABSTRACT 



The intermediate-mass star-forming core UYSO 1 has previously been found to exhibit intriguing features. While deeply embedded 
and previously only identified by means of its (sub-)millimeter emission, it drives two powerful, dynamically young, molecular 
outflows. Although the process of star formation has obviously started, the chemical composition is still pristine. We present Herschel 
PACS and SPIRE continuum data of this presumably very young region. The now complete coverage of the spectral energy peak 
allows us to precisely constrain the elevated temperature of 26-28 K for the main bulge of gas associated with UYSOl, which is 
located at the interface between the hot H n region Sh 2-297 and the cold dark nebula LDN 1657 A. Furthermore, the data identify 
cooler compact far-infrared sources of just a few solar masses, hidden in this neighbouring dark cloud. 
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1. Introduction 

Star formation occurs predominantly in structured environments 
in the case of intermediate and high-mass star formation. In 
particular, one often finds different stages of evolution located 
within close proximity. This poses a challenge for working on 
the earliest phases of star formation, which can only be revealed 
with long-wavelength observations. To observe these deeply em- 
bedded objects at the spectral peak of their emission, airborne 
or satellite missions are necessary, which in the past delivered 
only a very modest spatial resolution (>20")- We use ESA's 
new f ar-IR and sub-millimeter satellite Herschel (Pilbra tt et alj 
2010) — with its unprecedented spatial resolution — to scru- 
tinise a sample of very young low- and high-mass star-forming 
cores within our programme "Earliest Phases of Star Formation" 
(EPoS, PI: O. Krause). One of these targets is a core located in 
the Canis Majoris OB1/R1 region (e.g., Gregorio-Hetem et al. 
2009) of the outer Galaxy at a distance of ~ 1 kpc. Originally 
detected as a distin ct submillimetre sou rce in the vicinity of 
IRAS 07029-1215 (iForbrich et al.1 12004 . it was coyly named 
"UYSO 1" (unidentified young stellar object 1) since it was only 
detected at 450 and 850 //m and seemed to have no counterpart at 
other wavelengths. In the optical, the entire region is dominated 
by the strong Hn region Sh 2-297, excited by the early B star HD 
53623. The mid-IR emission is dominated by warm dust and 
PAH emission from the associated southern photon-dominated 
region (PDR). All these sources, however, do not coincide with 
UYSO 1, located more than 1' to the north-west. Further to the 



* Herschel is an ESA space observatory, with its science instru- 
ments PACS, SPIRE, and HIFI provided by European-led Principal 
Investigator consortia, and with important participation from NASA. 



west, one finds the optical dark cloud LDN 1657 A, about which 
little is known. Members of our group carried out a multiwave- 
length study of UYSO 1 (IForbrich etal.ll2009h . which revealed 
more puzzles. The source drives a dynamically young, but strong 
CO molecular outflow. The imprints of two crossed jets are vis- 
ible in shocked H2 emission in the near-IR, and their vertex is 
very close to the position of UYSO 1. Millimeter continuum 
interferometry resolved the source into two peaks separated by 
4'.'2. Spitzer/MlPS data indicated that at 24 and 70 //m, UYSO 1 
is still too deeply embedded to be directly detected, a finding we 
revise in the present paper for 70 fim with Herschel. In addition 
to the core being chemically pristine and showing no signs of 
more evolved chemistry, this object is very interesting for the 
investigation of the early evolution of (intermediate-mass) pro- 
tostars. 

2. Observations and data reduction 

The object UYSO 1 and its vicinity were observed in the blue (70 
//m), gree n (100 //m), and red (160/vm) filters with the PACS in- 
strument (Poglitsch et al. 2010) aboard Herschel on October 9, 
2009, within the science demonstration program. Scan maps in 
two orthogonal directions with scan leg lengths of 7' were ob- 
tained with a scan speed of 20"/s, for which the spatial resolution 
is 5'.'6, 6'.'8, and 1 1'.'3, respe ctively. The raw data were reduced 
with the HIPE software (e.g., Poglitsch et al.ll2010l) . version 3.0, 
build 455. Beside the standard steps leading to level- 1 calibrated 
data, a second-level deglitching and a correction for offsets in 
the detector sub-matrices were performed. Finally, the data were 
highpass-filtered, using a median window of the size of the full 
scan legs, to remove the effects of bolometer sensitivity drifts 
during data acquisition and 1/f noise. We masked out emission 
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Fig. 1. Left: Star-forming complex around UYSO 1 with PACS. Three-colour composite with 70 pm (blue), 100 /jm (green), and 
160 //m (red). Black cross at reference position: location of UYSO 1 derived by interferometry (see Fig. [2]i, white cross: peak of Sh 
2-297 cm emission, ellipse: error ellipse of IRAS 07029-1215 position. Right: SPIRE three-colour composite with 250 fim (blue), 
350 //m (green), and 500 //m (red). 



Table 1. Table of positions of relevant objects in the Herschel 
field. 



ID 


RA (J2000) 


Dec (12000) 




[ h : m : s] 


[":':"] 


IRAS 07029-1215 


07:05:16.9 


-12:20:02 


HD 53623 


07:05:16.75 


-12:19:34.5 


Sh-2 297 


07:05:15.4 


-12:19:39 


UYSO 1" 


07:05:11.14 


-12:19:00.2 


LDN 1657A-1* 


07:04:57.14 


-12:16:57.4 


LDN 1657A-2' 


07:04:58.03 


-12:16:50.8 


LDN 1657A-3* 


07:05:00.69 


-12:16:44.7 


LDN 1657A-4" 


07:05:03.50 


-12:16:33.7 


LDN 1657A-5* 


07:05:05.84 


-12:16:51.4 



"derived from astrometry-adjusted PACS 70 /jm map (see Sect. 13.1 
* derived from MIPS 24 pm map (see Sect.13^21 



Table 2. Point source fluxes (without colour correction). The 
overall uncertainties are 20% for all wavelengths except for 160 
pm (25%). 



ID 


F, [ly] 


F v 


[ly] 


F v [Jy] 


F v [Jy] 




(24 /im) 


(70 


pm) 


(100 pm) 


(160 pm) 


UYSO 1" 


- 




168 


239 


193 


LDN 1657A-1 


0.07 




0.6 


0.8 


1.5 


LDN 1657A-2 


0.10 




0.2 


0.2 


<0.2 


LDN 1657A-3 


0.09 




0.5 


0.9 


1.3 


LDN 1657A-4 


<0.002 




0.3 


1.1 


2.6 


LDN 1657A-5 


0.04 




0.4 


0.5 


0.7 



" Photometry obtained with large aperture (see Sect. 13. U . The corre- 
sponding SPIRE fluxes used for Fig. [3] are: (91.1±13.7) Jy (250 lira), 
(32.5±4.9) ly (350 pm), and (11.2+1.7) ly (500 pm). 



structures (visible in a first iteration) before computing and sub- 
tracting this running median; this step minimises oversubtraction 
of source emission. Finally, the data were projected onto a coor- 
dinate grid using th e photProject routin e inside HIPE, and flux 
calibration factors (Poglitsch et al. 2010) were applied. 

Maps at 250, 3 50, and 500 pm were obtained with SPIRE 
dGriffinetal .1120101) on October 19, 2009. Two 9' scan legs cov- 
ered the source with the nominal speed of 30"/s. The data were 



processed within HIPE with the standard photometer script up 
to level 1. During baseline removal, we masked out the high- 
emission area in the center. For these observations, no cross-scan 
data was acquired; therefore the resulting maps still show resid- 
ual stripi ng along the scan direction. An iterative destriping al- 
gorithm dBendo et al.ll2010l) was applied to mitigate this effect. 
We adopt FWHM beam sizes of 18'.'1, 2572, and 3679 at 250, 
350, and 500 pm, respectively. 



3. Results 

The Herschel data show the detailed structure of the entire star- 
forming complex containing UYSOl. The data are shown in 
Fig.Q] where we have combined the PACS and SPIRE data into 
two colour composites. In the southern region, there is a ridge of 
far-IR emission roughly oriented east-west. The extended emis- 
sion there covers the formal position of IRAS 07029-1215. The 
IRAS positions, however, have usually been derived from the 
IRAS 12 pm maps. We note that the reported position coin- 
cides with neither the much stronger far-IR peak associated with 
UYSO 1 (see below) nor the peak position of centimetre emis- 
sion arising from the Hn region (cf. FigJTJ. Parts of this emission 
ridge represent a PDR adjoining the Hn region. Nevertheless, the 
bright condensation east of the IRAS position within this PDR 
ridge might correspond to another embedded compact core. 

While this southern ridge is the brightest source at mid-IR 
wavelengths (A < 25 pm), the Herschel data present a far dif- 
ferent picture. The emission area around UYSO 1 dominates the 
Herschel bands: UYSO 1 itself is situated on top of a highly 
structu red emission plateau (not detected in the SCUBA 450 pm 
data of Forbri ch et al.ll200 4 j) with a similar ex tent and structure 
as the 850 pm map bv iForbrich et al.1 (12004]) . The emission in 
the close eastern vicinity (~ 30") of UYSO 1 is also a PDR ex- 
cited by the B star a s indicated by the [O i] line observed by 
IForbrich et al.l d2009l) . We do not delve into a detailed study of 
the te mperature distribution for the whole region (see Stutz et al. 
2010, for CB 244). However, from inspection of Fig. Q]it is ap- 
parent that the extended emission around UYSOl has its max- 
imum shortwards of 160 pm, while the region extending to the 
north-west contains emission that only becomes pronounced at 
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Fig. 2. PACS 70 yum image of the central region of UYSO 1 . 
The inset shows the PACS 70 /im PSF dPoglitsch etalJl2010h . 
rotated and scaled to match the UYSO 1 observation. The 
white contours denote the SMA 349 GHz continuum detection 
(H. Beuther, priv. comm.) 



the SPIRE wavelengths, indicating that even colder material is 
located in the LDN 1657 A dark cloud. 



3.1. UYS0 1 

At far-IR wavelengths, the UYSO 1 core dominates the contin- 
uum emission of the region. The PACS instrument offers un- 
precedented spatial resolution of around 5'.'6 for the 70 /mi filter 
providing the opportunity to relate th e far-IR emission to high- 
resolution data at other wavelengths. Forbrichet al. (2003) re- 
ported Spitzer/MlPS observations of UYSO 1 at 70 /mi as well as 
3 mm interferometric data pinpointing the two continuum peaks 
in UYSO 1 . A large position shift of > 9" was found between the 
mm interferometric sources and the 70 /im peak. Thus, the au- 
thors concluded that the optical depth at the column density peak 
is still too high for a direct detection of the sources at 70 /mi. 
This large discrepancy in peak positions is not confirmed by our 
PACS data. Figure [2] shows the PACS view on the central region 
around UYSO 1 . As a position reference for the two related mm 
sources, we use new SMA 349 GHz continuum data (H. Beuther, 
priv. comm.) that are in very good posi t ional agreement with the 
3 mm PdBI data from Forb rich et al.l d2009l) . The typical 1-cr 
absolute pointing accuracy of Herschel is 2"; however, devia- 
tions of up to 4" have been reported for a few operational days. 
Fortunately, the PACS data reveal several point sources north- 
west of UYSO 1 (not covered by the MIPS 70 pm data) that 
were also detected in the MIPS 24 pm map (see Sect. 13.21 ). We 
concatenated the PACS astrometry to the MIPS 24 pm astrom- 
etry (see Fig. [2j. As a result, the 70 pm peak is now close to 
the stronger of the two interferometric point sources. The two 
positions still do not perfectly coincide, but the remaining offset 
of ~ 2" is of the order of our final pointing uncertainties. Thus, 
there is a clear close association between the 70 /im peak and the 
millimetre source(s). As a test, we smoothed our PACS 70 /im 
data to the MIPS 70 /mi spatial resolution of 21" by convolving 
with an appropriate beam and regridding the resultant map to the 
same pixel scale. The final map is similar to the MIPS data, and 
the new peak position is indeed offset from that in the original 
PACS data by over 10". Because of the almost 4 times coarser 
spatial resolution of Spitzer at 70 //m, the true peak position had 
been smoothed out and merged with the surrounding PDR emis- 
sion, and thus appeared shifted. 

The comparison with the PACS 70 pm point-spread func- 
tion (PSF) in Fig. [2] illustrates that UYSO 1 is not an unresolved 
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Fig. 3. Single-temperature modified black-body SED fits left for 
UYSO 1, and right for LDN 1657A-4. Plotted on the same scale, 
this illustrates the different natures of the two regions. Both 
results were obtained by adopting the Ossenkopf & Henning 
( 1994) opacities (their Table 1 , /c-corumn 5). For T< 20 K, colour 
correction terms for PACS fluxes become larger than the mea- 
surement uncertainties (mainly at 70 pm). They were applied to 
the LDN 1657A-4 fluxes before the fitting. Plus signs: Herschel 
data, square: SCUBA (iForbrich et al.ll2004l) 
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Fig. 4. Compact sources in the LDN 1657A region (Table [TJ, 
north-west of UYSO 1. The cutouts show a region of 220" x 90" 
size at 24 pm (upper left), 70 /mi (upper right), 100 /im (lower 
left), and 160 //m (lower right). Note especially the steep SED 
rise of source 4, spanning from a non-detection at 24 pm to being 
the strongest source at 160 /im. 



point source. The PACS PSF has a characteristic triangular outer 
shape. However, UYSO 1 is marginally resolved at 70 //m; the 
slight elongation is along the connecting line between the two 
mm continuum sources. On the other hand, the stronger of the 
two outflo ws is also oriented in this north-west/south-east di- 
rection (cf. IForbrich et al.|l2009L their Fig. 2). Therefore, both 
source multiplicity and warm material along the outflow cone 
can explain our findings. 

The SED of UYSO 1 was estimated by fitting a Planck 
function, modified by wavelength-dependent opacities, to the 
Herschel data. Our approach was to use fluxes derived from 
uniformly sized apertures. We chose a circle with an area 
equivalent to the SPIRE 500 /im beam area (1543" 2 ). We 
used the opacities of coagul ated grains with thin ice mantles 
(Ossenkopf & Henning 1994). One temperature component was 
sufficient to obtain a very good fit (see Fig. [3]). The derived tem- 
perat ure of ~28.1 K was lower than the previous estimate of 40 
K in IForbrich eTail (l2009h . based on MIPS SED mode observa- 
tions that probably covered more of the adjacent warmer PDR 
material. Nevertheless, our photometry also includes a small 
contribution from the PDR, which is heated from the outside, 
thus, the intrinsic luminosity is slightly overes timated (in prin- 
cipal agreement with simulations performed bv lSteinacker et all 
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2006). Thanks to the clarified astrometry, actual Herschel fluxes 
for A < 450 yum could be included that exceed the Spitzer upper 
limits (which had been obta ined after s ubtracting an extended- 
emission component) of Forbrich et al. (2009). Consequently, 
the derived luminosity (244 L Q ) is higher than their estimate (50 
L Q ). We note that the derived total mass of 5.7 M strongly de- 
pends on the chosen dust model. We repeated the fitting with 
Draine & Lee (1984) opacities, deriving instead a total mass of 
26.4 M Q (and T = 26.1 K). 

3.2. New very cold source detections 

While UYSO 1 was the prime target of the Herschel observa- 
tions, we also discovered new sources 3' - 4' to the north-west 
of UYSO 1, within the neighbouring dark clou d LDN 1657 A . 
From the SCUBA 850 pm map shown in [Forbrich et al.l d2004l) . 
it was obvious that this region is associated with submillime- 
tre emission. The PACS data resolve this region into a "pearl 
necklace" of five compact sources (Table Q]). Nothing is known 
about these objects except the association of LDN 1657A-3 
with faint 2MASS K band emission. There is no counterpart 
at 1.4 GHz in the NVSS VLA survey. Unfortunately, the re- 
gion is at the very edge of previous ROSAT X-ray observations 
(Gregorio-Hetemet al. 2009) and thus not properly covered. 
However, inspection of the entire dataset of MIPS 24 pm dataj, 
of wh ich just the central part was published in iForbrich et al.l 
(2009), reveals that 4 of these 5 new sources are indeed detected 
(cf. Fig. [4j». At wavelengths <160yum, all these sources are un- 
resolved. Only at SPIRE wavelengths does a more pronounced 
extended emission component become apparent (Fig. Q]). Faint 
far-IR filaments seem to connect UYSO 1 and these northwest- 
ern sources. Hence, we assume that they are associated with the 
general star-formation complex. The distances between neigh- 
bouring LDN 1657A sources are quite similar (39" -43"), only 
the two western sources 1 and 2 being closer together (15"). 
At a distance of 1 kpc, this translates into around 0.2 pc (and 
0.073 pc, respectively) of projected distance. The distance 0.2 
pc is roughly the Jeans length of a cold (16-18 K), medium- 
dense (1 -2xl0 4 crrT 3 ) gas. Thus, these sources most likely rep- 
resent distinct cores. We report the P ACS photometry in T abled 
based on PSF photometry (see also Henni ng et al.ll2.Q10h using 
the Starfinder tool (Diolaiti et al. 2000) and reference PSFs pro- 
vided by the PACS instrument team. Most of these sources are 
detected at 24 yum, which is indicative of a second, warmer dust 
component. 

Source LDN 1657A-4 has the steepest rising SED slope 
of all these sources (Figs. [3j |4|. At 24 pm, the point source 
was undetected; we instead found a 24 yum shadow at its loca- 
tion, caused by the high column density of the obscuring ma- 
terial associated with t he cold core. We used the approach of 
IVasvunina et alj (120091) to transform the 24 //m shadow con- 
trast into a column density N(H). The intensity ratio of the 
surrounding extended emission to the darkest pixel in the 24 
/mi shadow was 1.51, after subtraction of the estimated zodi- 
acal foreground of 25.56 MJy/sr. The optical depth was then 
r = In 1.51, and the column density was given by N(H) = t/ct. 
We assumed t he extinction cross-s ection <x from the Ry-5.5B 
dust modejjof Weingartner & Draine (2001), appropriate for the 
denser shielded regions of molecular clouds. We estimated a col- 
umn density of 3.2 x 10 22 crrT 2 . This corresponds to Ay -26 mag, 



using the lRvterl (1 19961) calibration. A more el aborate modellin 
of the fore- and background contributions (e.g. JStutz et 



lodelling 
al.l2007l) 



1 The region has not yet been observed with Spitzer at wavelengths 
< 24 yum. The compact sources are not detected in the MSX images. 

2 1 .29x 1CT 23 cnr/H for MIPS 24 isophotal wavelength of 23.68 pm. 



would probably increase the calculated column density, but is 
beyond the scope of this letter. In addition, the standard (coarse) 
pixel scale of 2'.'45 in the MIPS 24 pm data used as well as an 
increased level of striping in the data limit the shadow contrast. 
Thus, the derived 24 pm shadow column density represents a 
lower limit. Using the SPIRE 250 pm peak flux measured at the 
position of LDN 1657A-4 (7 23 Jy/ 18"1 beam), the values of 
Ossenkopf & Henning (1994) opacities at this wavelength (18 
cm 2 /g), and the temperature derived from the SED fit to the pho- 
tometry data (17 K, see Fig.[3j right), we also calculate a column 
density of 3.2 x 10 22 crcT 2 . This value is certainly a lower limit, 
since the relatively large SPIRE beam smears out the true col- 
umn density peak. Nevertheless, the very good agreement be- 
tween the two independent methods is compelling. Finally, in- 
tegration of the 24 pm extinction map over the SPIRE 250 pm 
beam area results in a total mass of 1.6 M , a higher value than 
that derived from SED fitting (0.8 M ) but still comparable. 

4. Conclusions 

We have presented new Herschel PACS and SPIRE scan map 
data for the star-forming complex containing the intermediate- 
mass core UYSO 1 . They show the dust emission structures as- 
sociated with this region in unprecedented detail. The high spa- 
tial resolution of the PACS 70 pm map facilitates the differen- 
tiation of the true emission peak position from emission aris- 
ing from the surrounding PDR material. The PACS data show 
that the 70 pm emission is closely associated with the known 
location of the central millimetre peak(s) in UYSO 1 . This re- 
vises an earlier finding of large offsets (based on Spitzer/MlPS 
70 pm data). Hence, the measured Herschel photometry can be 
used with more confidence than for previous data to estimate 
the SED of UYSO 1 . Since the core was not clearly detected at 
24 pm (which otherwise would indicate a distinct component of 
warmer dust), a one-temperature modified black-body compo- 
nent with T = 28 K fits the data between 70-500yum very well. 
While clearly not in the hot-core/hot-corino regime, this temper- 
ature is higher than the typical values of < 1 7-20 K found in other 
very young regions of star formation. The known outflow activ- 
ity can be one explanation. However, t he influence of the n earby 
PDR has also to be taken into account ( Forbrich et al. 2009). The 
Herschel data reveal the serendipitous discovery of five compact 
very red sources, situated to the north-west of UYSO 1, within 
the neighbouring dark cloud LDN 1657 A. These objects have 
lower temperatures than UYSO 1 . A particularly intriguing ob- 
ject is LDN 1657A-4, which has a very red SED and exhibits 
a 24 pm shadow. We note that two different methods for the 
column density e stimation provided very consis tent results. For 
both objects, the lOssenkopf & Henning (1 19941) grain opacities 
seem more suitable than the use of ISM grains. 
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